Role of non-collective excitations in low-energy heavy-ion reactions 
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We investigate the effect of single-particle excitations on heavy-ion reactions at energies near 
the Coulomb barrier. To this end, we describe single-particle degrees of freedom with the random 
matrix theory and solve the coupled-channels equations for one-dimensional systems. We find that 
the single-particle excitations hinder the penetrability at energies above the barrier, leading to a 
smeared barrier distribution. This indicates that the single-particle excitations provide a promising 
way to explain the difference in a quasi-elastic barrier distribution recently observed in 20 Ne + 
90 ' 92 Zr systems. 
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I. INTRODUCTION 

Heavy-ion reactions near the Coulomb barrier provide 
a good opportunity to investigate an interplay between 
the reaction process and internal excitations in the col- 
liding nuclei. For example, it is well known that sub- 
barrier fusion cross sections are significantly enhanced as 
compared to a prediction of a simple barrier penetration 
model because of the couplings of the relative motion be- 
tween the colliding nuclei to nuclear intrinsic degrees of 
freedom [l| . It has been well recognized by now that the 
enhancement of fusion cross sections can be explained in 
terms of a distribution of Coulomb barrier heights orig- 
inated from the couplings [2|. The barrier distribution 
can be actually extracted directly from the experimental 
data by taking the second derivative of the product of 
fusion cross section <7f us and the center of mass energy E 
with respect to E, that is, d 2 (Ea {vLS (E)) /dE 2 @j. The 
experimental data have clearly shown that the barrier 
distribution makes a useful representation to understand 
the reaction dynamics of heavy-ion subbarrier fusion re- 
actions ply]- 

A similar concept of barrier distribution has been ap- 
plied also to heavy-ion quasi-elastic scattering (that is, 
a sum of elastic, inelastic scattering, and transfer reac- 
tions) [1, Q. In this case, the barrier distribution is de- 
fined as the first derivative of the ratio of the quasi-elastic 
scattering cross section at a backward angle <7 qG i to the 
Rutherford cross section or, that is, —d(a qe \/aB,)/dE. 
The fusion and quasi-elastic barrier distributions have 
been found to behave similarly to each other at least in 
a qualitative way [J, Q . 

In order to analyse the subbarrier enhancement of fu- 
sion cross sections and fusion and quasi-elastic barrier 
distributions, the coupled-channels method has been suc- 
cessfully employed. Typically, a few low-lying collective 
excitations are taken into account in a calculation |l|, |7j. 
With this approach, the barrier distribution arises natu- 
rally through the eigen-channel representation [J, fsl— FlOJ| . 

Recently, however, a few experimental data which 
cannot be accounted for by the conventional coupled- 
channels calculation have been obtained [FU - Rol ]. One of 



these examples is the quasi-elastic scattering experiment 
for the 20 Ne+ 90,92 Zr systems |15| ■ The experimental data 
show that the quasi-elastic barrier distribution for these 
systems behaves in a significantly different way from each 
other: the barrier distribution for the 20 Ne+ 92 Zr sys- 
tem is much more smeared than that for the 20 Ne+ 90 Zr 
system [15|]. On the other hand, the coupled-channels 
calculations that take into account the collective rota- 
tional excitations in 20 Ne as well as the vibrational exci- 
tations in 90 > 92 Zr lead to similar barrier distributions for 
both systems, because the strongly deformed 20 Ne nu- 
cleus mainly determines the barrier structure while the 
difference in the collective excitations in the two Zr tar- 
gets plays a minor role. In Ref. [15|, it was suggested 
that single-particle excitations in the colliding nuclei are 
responsible for smearing the barrier distribution for the 
20 Ne + 92 Zr system. Notice that the single-particle ex- 
citations are expected to be much more important for 
the 92 Zr nucleus compared to the 90 Zr nucleus, which 
has the N — 50 shell closure. In fact, while there are 
only 12 states in the 90 Zr nucleus up to 4 MeV, there 
are 53 known states in the 92 Zr nucleus [lg. For 5 MeV, 
the number of known states is 35 and 87 for 90 Zr and 
92 Zr, respectively (for another nucleus, 116 Sn, there are 
81 known levels up to 3.9 MeV and 112 levels up to 4.3 

MeV [IEGJ]). 

The aim of this paper is to investigate the effect of low- 
lying single-particle excitations on low-energy heavy-ion 
reactions, as conjectured in Ref. 15[. In order to un- 
derstand qualitatively the effect of non-collective excita- 
tions, in this paper we shall use a schematic model, that 
is, one-dimensional barrier penetration in the presence of 
the couplings to intrinsic degrees of freedom. The single- 
particle degrees of freedom can be described in several 
ways [13-[2J|. For instance, Ref.|19| used the Lindblad 
approach to discuss the role of quantum decoherence in 
deep subbarier hinderance of fusion cross sections. In this 
paper, we employ the random matrix theory (RMT) to 
describe the single-particle excitations (see Refs. [25 
|28| for recent reviews on RMT). The random matrix 
approach for heavy-ion reactions has been developed in 
1970's by Weidenmiiller and his collaborators in order 
to analyse heavy-ion deep inelastic collisions (DIC) [20l — 



|23[ . At that time, they derived the transport coefficients 
based on RMT [29j and solved the classical transport 
equations (see also Refs. [30,[3l|). The RMT has also 
been employed to discuss quantum dissipation [32h34| . In 
this paper, instead of solving the classical equations, we 
directly solve the coupled-channels equations quantum 
mechanically by including the single-particle excitations 
described by RMT. Our approach is therefore similar to 
that in Ref. [24| , in which the coupled-channels equations 
with 200 dimension were solved for a one-dimensional 
model us ing a semi-classical approximation. In contrast 
to Ref. [24| . we apply our formalism to the subbarrier 
regime without using the semi-classical approximation. 
This will enable us to assess the effect of single-particle 
excitations on quantum tunneling and thus on the bar- 
rier distribution. By treating the single-particle states 
explicitly, we can also discuss the excitation spectra as a 
function of incident energy. 

The paper is organized as follows. In Sec. II, we de- 
tail the coupled-channels formalism with single-particle 
excitations described by RMT. In Sec. Ill, we apply the 
formalism to one-dimensional models for quantum tun- 
neling. We discuss the effect of single-particle excita- 
tions on the barrier penetrability, the barrier distribu- 
tion, and the excitation spectra. Using the results for 
the one-dimensional model, we also discuss the effect of 
non-collective excitations on quasi-elastic barrier distri- 
bution for the 20 Ne+ 92 Zr system. We then summarize 
the paper in Sec. IV. 



II. FORMALISM 

A. Coupled-channels method 

The aim of this paper is to discuss the effect of single- 
particle excitations on one-dimensional barrier penetra- 
bility. For this purpose, we assume the following Hamil- 
tonian: 

H = -^^2+ y rcl(20 + #0(0 + Koup(x, 0- (1) 

Here, fi is the reduced mass and V le i(x) is a potential for 
the relative motion. Hq(£) is a Hamiltonian for the in- 
trinsic degrees of freedom of the colliding nuclei, and the 
last term, V coup (x,t;), is a coupling Hamiltonian between 
the relative motion and the internal degrees of freedom. 
The coupled-channels equations for this Hamiltonian 
are obtained by expanding the total wave function in 
terms of the eigen functions of Hq(^) and read, 



H 2 d 2 



2/i dx 2 



+ V re \(x) + e n - E \ 4> n (x) 

+ ^2v nm (x)^ m {x) = 0. (2) 



Here, e„ and ip n {x) are the excitation energy and the 
wave function for the n-th channel, respectively. V nm (x) 



is a coupling matrix and is a function of the coordinate 
x. 

The coupled-channels equations are solved by imposing 
the boundary conditions of 

ip n (x) -> <5„, e~ tkoX + r„ e iknX for x -> +oo (3) 
-> t n e ~ lknX for x -> -00, (4) 



where k n = y 2/x(2? — e n )/h 2 is the wave number for 
the n-th channel, and represents the entrance channel. 
We have assumed that the projectile is incident from the 
right hand side of the potential barrier. With the trans- 
mission coefficients t n , the penetration probability for the 
inclusive process is calculated as 



"ii 



p(e) = J2 p ^ e ) = Et^ 



(5) 



The barrier distribution is obtained by taking the deriva- 
tive of P(E), that is, dP{E)/dE % 

In order to take into account the single-particle exci- 
tations, as we will show in the next section, one has to 
include a large number of channels. Since it is time and 
memory consuming to solve the coupled-channels equa- 
tions with a large dimensionality, in this paper we employ 
a constant coupling approximation [8J. In this approxi- 
mation, the coupling matrix is assumed to be a constant 
over the whole range of x. Then, one can diagonalize the 
matrix A = (V nm + £ n 3n,m) with a coordinate indepen- 
dent unitary matrix U, 



UAU^ =diag{X 1 ,X 2 ,---} 



(0) 



where, Ai, A2, • • • are the eigenvalues of A. Transforming 
the channel wave functions as 



^n(x) 






U nmtym \% ) 3 



(7) 



the coupled-channels equations are transformed to a set 
of the uncoupled equations, 

h 2 d 2 1 ~ 

- — -^+V iel {x) + \ n -E\ip n (x)=0- (8) 

We call the transformed channels the eigen-channels and, 
for each eigen-channel, V n (x) = V,. c i(a;) + A n the eigen- 
potential. 

The boundary conditions for ip n (x) are given by 

ij> n {x) -» U n0 ( e -' lkox + 7 n e lkox ) for x -> +00 (9) 

and 

Mx) -> U n0 t n e~ lkoX for x -> -00 (fO) 

where the reflection and transmission coefficients are re- 
lated to the original coefficients in Eqs. (j3]) and (j4]) by 



\ U nm r m /U nQ and t„ 






U nm t m /U nQ , respec- 



tively. Here, we have assumed that the excitation ener- 
gies are small compared to the incident energy so that 



k n can be approximated by ko [8[ . Using the coefficients 
t n , the penetrability is calculated as 



P(E)=Y,\ U -o\ 2 \tn(E)\' 



(11) 



The reflection coefficients in the original basis are given 
by 



22(U l )nmU m0 ''i 



(12) 



from which the Q-value distribution (that is, the excita- 
tion spectrum) is computed as 

/w = ErM 2<5 ( e - e ")~EM 2 <^- £ «)- ( 13 ) 



B. Coupling matrix elements 

We solve the coupled-channels equations, (2), using the 
constant coupling approximation by including both col- 
lective and non-collective excitations. For the collective 
excitations, we assume either the vibrational or the rota- 
tional couplings. The coupling matrix for the vibrational 
coupling is given by 



(V nm ) 



1 

1 



(14) 



if we truncate the phonon space up to 1-phonon state [7| . 
Here, F is a coupling constant, and we have assumed the 
linear coupling. For the rotational coupling, the coupling 
matrix is given by 



(V n 




2\/5 


6 


7 


7 


(i 


20\/5 


7 


77 - 





1 


6 


20\/5 


7 


77 


20v^ 


486 


77 


1001 



(15) 



up to the 4 + state in the rotational bandpj, where i<2 
and Fi are the quadrupole and hexadecapole coupling 
strengths, respectively. 

For the single-particle excitations, we consider an en- 
semble of coupling matrix elements based on the random 
matrix theory [20l - l22l |. We assume that the matrix ele- 
ments are uncorrelated random numbers obeying a Gaus- 
sian distribution with zero mean. That is, we require that 
the first and the second moments of the coupling matrix 
elements satisfy the following equations |35j 



V nm {x) = 



WQ _(*n-*m) 2 



Qnm — 



v /p(e„)p(e m ) 



(16) 
(17) 

(18) 



where the overbar denotes an ensemble average and p(e) 
is the nuclear level density. Here, we have assumed the 
coordinate independent matrix elements according to the 
constant coupling approximation. 

For the single-particle excitations, we generate the 
coupling matrix elements according to these equations 
many times. For each coupling matrix, we do not vary 
the matrix elements for the collective excitations, which 
are uniquely determined once the coupling is specified. 
For each coupling matrix, we solve the coupled-channels 
equations and calculate the penetrability and the reflec- 
tion probability. The physical results are then obtained 
by taking an average of these quantities. 

In the actual calculations shown in the next section, 
we discretize the quasi-continuum single-particle spec- 
trum in the coupled-channels equations [36|, |37( (see also 
Ref.|38j). Introducing the level density p(e), 



p( e ) =X^( e-e ™)' 



(19) 



the coupled-channels equation can be written in the fol- 
lowing form : 



h 2 d 2 

-75- -J-"2 + Viel( x ) + e n~ E \ 4>n{x) 



+ / dep(e)V ne (x)ip e (x) = 0. (20) 

In this equation, we assume a quasi-continuum spectrum 
for the single-particle excited states, and discretize the 
integral with a constant energy spacing, Ae. For the 
ground state and the collective excitation channels, we 
then obtain, 



- E V nm(x)lp m (x) 

m^sp 

- J2 Ae P(£rn)V nem (x)lP em (x) = 0, (21) 
m(Esp 



while for the single-particle channels denoted by e n we 
obtain 



h 2 d 2 1 

~2"fa? + Vicl ^ + €n ~ E I ^ £ "^ 

- 2_j V ^n(x)lp m (x) 

m^sp 

- J2 Aep(e m )V enem (x)il> em (x) - 0. (22) 

mEsp 



Here, m <£ sp denotes a summation over the ground state 
and the collective channels, while m £ sp a summa- 
tion over the single-particle channels. These equations 



can b e express ed in a simpler way by multiplying a fac- 
tor y/p(e n )Ae for each index e n representing the single- 
particle channels. That is, 



III. RESULTS 



A. Vibrational coupling 



^ e „ (x) = y/p(e n )Ae ^„ (x), (23) 

V n£m (x) = ^Jp{e m )Ae V nem (x), (24) 

V enem (x) = Aey/ p(e n )p(e m ) K„ em (a). (25) 

With these wave functions and the coupling matrix ele- 
ments, Eqs. <J2U) and J52J read, 



+ V TO i(x) + e n - E \ ip n (x) 



2(i dx 2 

+ Y, V nm^m(x) + J^ ^ m (x)A m {x) = 0, (26) 
m^sp 



mtEsp 



and 



+ V^{x)+e n -E\^ tn {x) 



2p dx 2 

+ E V^rn(x) + J2 Ve n e m {x)^e m {x) = 0, (27) 
m£sp roSsp 



respectively. From Eqs. ([3]) and (U), the boundary con- 
ditions for il) €n (x) are given by 



tpe n (x) -> y/p{e n )Ae t n e lk 

-ik n x r 



CwiC 



r ,x — > — oo 



(28) 



and 



^e„(a;) -> y/p(e n )Ae r n e l 
= r n e lknX for x 



(29) 



where, £„ = N /p(e„)Ae £„ and r„ = ^p(e„)Ae r„. The 
penetrability is then given by 



¥) = ErW 2+ Er^) A£ l' 



n^sp 

E 



feo 



n£sp 



fc 



2 

n\ -i 



\tn\ 2 + 



E 

nGsp 



fco 



(30) 



We now numerically solve the coupled-channels equa- 
tions and discuss the effect of single-particle excitations 
on barrier penetrability. We first consider the vibrational 
coupling. We assume that there is a collective vibra- 
tional state at 1 MeV, whose coupling to the ground 
state is given by Eq. ([11]) with F = 2 MeV. For 
the single-particle states, we consider a level density 
given by p(e) = p e 2 ^ 1 with p = 0.039 MeV -1 and 
a = 29/8 MeV -1 , starting from 2 MeV. The value of 
Po was determined so that the number of single-particle 
levels is 200 up to 5 MeV. For the parameters for the 
couplings in Eq. (fT8|) . we follow Ref. [35| to use A=7 
MeV. We arbitrarily choose the coupling strength to be 
u)q = 0.005 MeV. The energy spectrum for this model is 
shown in FigfTJ For the potential for the relative motion, 
V Te i(x), we use a Gaussian function 



V rel (x)=V B e *%. 



(31) 



with Vb = 100 MeV and Sq = 3 fm Q. The reduced mass 
p is taken to be 29mjv, m^r being the nuclcon mass. 
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FIG. 1: (Color online) The energy spectrum for the model 
calculation we employ. There is a collective vibrational state 
at 1 MeV, while single-particle states exist from 2 MeV with 
an exponentially increasing level density. 



This method with a constant energy step considerably 
reduces the computation time as compared to the case of 
treating the exponentially increasing number of single- 
particle levels as they are. It also validates the use of 
the random matrix theory. An important assumption 
in RMT is that the ensemble average of a quantity is 
equivalent to the energy average of that quantity over 
the spectrum [26j. In our case, we expect that the en- 
semble average of the calculated results corresponds to 
the energy average of the same quantities within the en- 
ergy spacing of Ae. 



Figure [2] shows the penetrabilities thus obtained. The 
corresponding barrier distributions are shown in Fig. [3] 
The dotted and the dashed lines show the results without 
the channel couplings and those only with the collective 
excitation, respectively. The solid line shows the results 
with both the single-particle excitations and the collec- 
tive excitation. We include the single-particle states up 
to e max =23 MeV with energy spacing of Ae=0.02 MeV. 
With this model space, the number of channels included 
is 1013 (we treat the low-lying single-particle states as 



discrete states when the energy spacing is larger than 
Ae). This result is obtained by generating the coupling 
matrix elements 30 times to take an ensemble average. 
We have found that the fluctuation around the average is 
small. For instance, at E = 1 00 MeV, the averaged pene- 
trability is P — 0.622, whereas the root-mean-square de- 
viation is 4.366xl0~ 3 . As we mentioned in the previous 
section, the results shown in Fig. 2 are obtained with the 
constant coupling approximation. With a smaller value 
of e max , we have solved the coupled-channels equations 
exactly and have confirmed that the constant coupling 
approximation works qualitatively well. 

The collective excitation leads to a double peaked 
structure of barrier distribution. One can see that the 
single-particle excitations suppress the penetrability at 
energies above the barrier, and at the same time smear 
the higher energy peak in the barrier distribution, al- 
though the main structure of the barrier distribution is 
still determined by the collective excitation. The single- 
particle excitations also lower the barrier and thus in- 
crease the penetrability at energies below the barrier, due 
to the well-known potential renormalization [391 ]. 

The Q- value distribution for the reflected flux is shown 
in Fig. [4] at four incident energies indicated in the figure. 
For a presentation purpose, we fold the discrete distribu- 
tion with a Lorentz function, 



K £ ) = z 



V 



>1 



(32) 



with the width of r\ =0.2 MeV. That is, with the function 
defined by Eq. (fl"3|) . we compute 



(33) 



F(E*) = J def(e)g(E* - e) 

= £1 



7T (E* - e n ) 2 + rf ' 



In the figure, the peaks at E* = MeV and E* = I MeV 
correspond to the elastic channel and the collective exci- 
tation channel, respectively. One can see that at energies 
well below the barrier the elastic and the collective peaks 
dominate in the distribution. As the energy increases, the 
single-particle excitations become more and more impor- 
tant. This behaviour is consistent with the experimental 
Q-value distribution observed for 16 O+ 208 Pb 0,IH1 and 
16 0+ 184 W [41| reactions. At energies above the barrier, 
the single-particle contribution is even larger than the 
contribution of the elastic and the collective peaks. 



B. Rotational coupling 

Let us next consider the rotational coupling. For this 
purpose, we mock up the 20 Ne + 92 Zr system. That is, 
we consider the rotational excitations in 20 Ne up to the 
4 + state and the single-particle excitations in 92 Zr. The 
energies of the rotational states are thus e 2 + = 1-634 
MeV and e 4 + = 4.248 MeV for the 2+ and 4+ states, 
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FIG. 2: (Color online) The potential penetrability 
obtained with several methods. The dotted line is 
obtained without channel coupling, while the dashed 
line takes into account only the collective vibrational 
excitation. The solid line shows the result with both 
the collective and the single-particle excitations. 
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FIG. 3: (Color online) The barrier distribution de- 
fined by the first derivative of the penetrability. The 
meaning of each line is the same as in Fig[2] 



respectively. The values of the coupling strengths F% and 
Fa in Eq. (|15|) are estimated with the collective model for 
the coupling form factor at the barrier position with the 
deformation parameters of /?2 = 0.46 and f}± = 0.27. This 
yields F 2 = -6.892 MeV and F 4 = -4.632 MeV. For the 
single-particle excitations in 92 Zr, we consider the energy 
range of 2 MeV < e < 16 MeV, with the exponential level 
density with p = 0.039 MeV^ 1 and a = 30/8 MeV" 1 . 
For the coupling strength, we use wq = 0.005 MeV and 
A = 4.0 MeV. These parameters are adjusted so that the 
rotational excitation in 20 Ne gives the main structure of 
the barrier distribution. 

In the calculation shown below, we also include the 
mutual excitations of the projectile and the target nuclei. 
In order to avoid closed channels, we introduce the energy 
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FIG. 4: (Color online) The Q-value distribution for the re- 
flected flux at four energies as indicated in the figure. It is ob- 
tained by smearing the discrete distribution with a Lorentzian 
function with the width of 0.2 MeV. The peaks at E*=0 MeV 
and 1 MeV correspond to the elastic and the collective exci- 
tation channels, respectively. 



cutoff and include those channels whose total excitation 
energy is below 16 MeV. With this set-up, the number of 
channels included is 1688. 

We use the Gaussian function for the potential for the 
relative motion with Vb = 51.76 MeV and so = 2.475 fm. 
This yields the same barrier height and the curvature as 
those with a Woods-Saxon potential with the parameters 
of V = 59.9 MeV, r = 1.2 fm, and a = 0.63 fm for the 
20 Ne + 92 Zr system. The reduced mass is taken to be 
20x 92 tojv/112. 

Figs. [5] and |5] show the penetrability and the barrier 
distribution for this model, respectively. The meaning of 
each line is the same as in Figs. [2]and[3]for the vibrational 
coupling. With only the collective rotational excitations, 
there are three eigenbarriers whose height is 48.2, 53.3, 
and 55.6 MeV. In contrast to the vibrational coupling, for 
the rotational coupling with a prolate deformation the 
main peak in the barrier distribution is not the lowest 
energy peak. For the parameters we use, the highest 
energy barrier carries a relatively small weight and the 
barrier distribution has only two visible peaks. The effect 
of single-particle excitations on the barrier distribution is 
similar to the case for the vibrational coupling and smears 
the higher energy peak in the barrier distribution. 
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FIG. 5: (Color online) The same as Fig. [2] but for 
the rotational coupling. 
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FIG. 6: (Color online) The same as Fig. [3] but for 
the rotational coupling. 



C. Quasi-elastic barrier distribution 

Using the eigenbarriers and the corresponding weight 
factors obtained in the previous subsection, one can com- 
pute the quasi-elastic scattering cross sections and the 
quasi-elastic barrier distribution in a three-dimensional 
space. That is, in the eigen-channel representation, the 
quasi-elastic scattering cross section is given by [a, [42( , 



<J qe i(E 7 9) =22 w l <J c \{E - X t 



(34) 



where wi — \Uio\ 2 is the weight factor for the i-th eigen- 
channel and <7 e i is the elastic scattering cross section. In 
order to calculate the elastic scattering cross sections, we 
use the Woods-Saxon potential indicated in the previous 
subsection. For the imaginary part of the optical po- 
tential, we assume an internal absorption, in which the 
imaginary part is well localized only inside the Coulomb 
barrier. 
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FIG. 7: (Color online) The quasi-elastic barrier distribution 
for the 20 Ne+ 92 Zr system calculated with the weight factors 
and the eigenbarrier heights for the one-dimensional rota- 
tional model discussed in Sec. Ill B. The meaning of each 
line is the same as in Fig. [3] 



Fig. [7] shows the quasi-elastic barrier distribution ob- 
tained with the same Wi and A» as in the previous sub- 
section for the one-dimensional model. We use the point 
difference formula with A25 c . m .= 2 MeV to calculate the 
quasi-elastic barrier distribution. The meaning of each 
line is the same as in Fig. [3] Because the quasi-elastic 
barrier distribution is by itself smeared more than the fu- 
sion barrier distribution 5] , and also because we use the 
point difference formula rather than taking the deriva- 
tive, the higher energy peak in the barrier distribution is 
more smeared by the single-particle excitations as com- 
pared to the one-dimensional calculation shown in the 
previous subsection. The difference between the dashed 
line (the collective excitations only) and the solid line 
(the collective -I- single-particle excitations) is similar 
to the difference in the experimental quasi-elastic bar- 
rier distribution between 20 Ne+ 90 Zr and 20 Ne+ 92 Zr sys- 
tems. We therefore conclude that the single-particle ex- 
citations indeed provide a promising way to explain the 
difference in the quasi-elastic barrier distribution for the 

20 Ne+ 90,92 Zr systems _ 



applied the model to one dimensional barrier penetration 
problems by using the constant coupling approximation. 
In addition to the single-particle excitations, we also in- 
cluded the collective excitations with cither a vibrational 
or a rotational character. We calculated the potential 
penetrability, the barrier distribution and the Q-value 
distribution. We also calculated the quasi-elastic bar- 
rier distribution for the 20 Ne + 92 Zr system using the 
eigenbarriers and their weight factors obtained with the 
one-dimensional model. Our calculations show that the 
single-particle excitations hinder the penetrability at en- 
ergies above the barrier and smear the high energy part 
of the barrier distribution. In the Q-value distribution, 
we found that the contribution from the single-particle 
excitations increases significantly as the incident energy 
increases. 

The experimental quasi-elastic barrier distributions are 
considerably different between 20 Ne + 90 > 92 Zr systems, 
despite that the coupled-channels calculations with col- 
lective excitations in the colliding nuclei lead to similar 
barrier distributions to each other. Our calculations im- 
ply that the difference can be indeed accounted for by 
the non-collective excitations in the target nuclei, as has 
been conjectured in Ref. [15| 

In order to make a quantitative comparison to the ex- 
perimental data and draw a definite conclusion on the 
quasi-elastic barrier distribution for the 20 Ne + 90 > 92 Zr 
systems, it will be an interesting future work to extend 
our study presented in this paper to three-dimensional 
calculations without resorting to the constant coupling 
approximation. This is so espeically because the constant 
coupling approximation which we employed in this paper 
may introduce a significant phase error and thus leads to 
an inconsistent angular dependent interference between 
different partial waves. It may also be interesting to see 
whether this model accounts for the hindrance of fusion 
cross sections at deep sub-barrier energies recently found 
in several systems. For this purpose also, we would have 
to take into consideration the coordinate dependence of 
the coupling form factor, especially around the touching 
point of the colliding nuclei |43|. 
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